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ABSTRACT 


Arrival-time  estimation  by  adaptive  thresholding  is  described. 
The  probability  density  of  ai  rival  time  is  derived  for  differentiable 
Markov  processes.  The  special  case  of  additive,  stationary  noise  is 
given  particular  attention.  A  direct  derivation  of  the  probability 
density  of  arrival  time  for  pulses  with  sharply  rising  edges  is  given  for 
arbitrary  noise.  1116  results  are  applied  to  the  Gaussian  and  Rice 
distributions.  Comparison  with  the  Cramer-Rao  bound  of  estimation 
theory  indicates  the  asymptotic  optimality  of  adaptive  thresholding 
for  these  two  distributions. 
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ARRIVAL- TIME  ESTIMATION  BY  ADAPTIVE  THRESHOLDING 


INTRODUCTION 

One  me+nod  of  determining  the  arrival  time  of  a  signal  is  to  measure  the  time  at  which 
the  signal  crosses  a  certain  threshold  level.  Figure  1  shows  a  puLe  with  arrival  time  defined 
to  be  the  time  at  which  the  pulse  exceeds  level  A .  Suppose  a  series  of  pulses  are  received. 
In  a  practical  system  the  amplitude  E  will  usually  vary  from  pulse  to  pulse.  If  the  rise 
time  does  not  change,  the  arrival  time  relative  to  the  leading  edge  will  vary  from  pulse  to 
pulse,  even  if  no  noise  is  present.  To  remedy  the  situation,  adaptive  thresholding  can  be 
used.  Adaptive  thresholding  is  defined,  to  be  arival-time  estimation  by  thresholding  in 
which  the  level  A  is  always  a  fixed  fraction  of  E.  Typically,  A/E  is  set  equal  to  1/2. 


*«  TIME 

Fig.  i—Received  pulse  with  a  threshold  crossing 


A  practical  adaptive-thresholding  system  is  shown  in  Fig.  2.  The  received  signal  is 
split  into  two  branches.  In  the  lower  branch  the  signal  is  first  delayed  by  an  amount  ex* 
ceeding  its  rise  time,  which  is  assumed  to  be  known.  Then  the  delayed  signal  passes 
through  an  amplifier  wit!  a  gain  of  G.  The  latter  is  a  negative  quantity  such  that  iCj"5  . 
equal  to  the  desired  val'  .e  of  A  'E.  The  amplifier  output  is  added  to  the  original  received 
s'gnal.  The  adder  output  crosses  the  zero  level  at  a  time  equal  to  the  time  at  which  the 
level  A  is  crossed  plus  the  constant  delay.  In  the  important  special  case  of  a  half-ampli¬ 
tude  arrival  time,  A/E  -  1/2  ?md  G  =  -2.  A  detailed  analysis  of  this  adaptive  thresholder 
and  an  alternative  one  can  be  found  in  the  literature  (1). 

Arrival-time  estimation  by  adaptive  thresholding  is  not  the  optimum  estimation  pro¬ 
cedure  (2-6).  However,  it  -will  be  shown  that  adaptive  thresholding  provides  an  asymp¬ 
totically  efficient  estimate  at  least  in  some  important  cases.  Furthermore,  the  simplicity 
of  implementing  adaptive  thresholding  systems  gives  the  method  practical  importance. 
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An  approximate  treatment  of  adaptive  thresholding  in  the  presence  of  Gaussian  noise 
lias  previourly  been  done  (1).  In  this  report  the  general  case  of  differentiable  Markov 
processes  is  studied. 


GENERAL  FORMULATION 

Figure  3  shows  an  exaggerated  example  of  a  multiple  threshold  crossing,  which  may 
occur  when  noise  corrupts  die  received  signal.  The  combined  signal  and  noise  is  denoted 
by  r(f).  When  multiple  crossings  occur  dose  together,  there  is  some  ambiguity  as  to  which 
crossing  represents  the  actual  pulse  position,  i.e.,  which  crossing  is  closest  to  that  crossing 
which  exists  in  tire  absence  of  noise.  A  somewhat  arbitrary  but  simple  strategy  is  for  the 
adaptive  thresholder  to  detect  only  the  first  level  crossing  for  each  pulse  it  receives.  Thus 
we  seek  to  derive  the  probability  density  of  the  first  threshold  crossing  time  of  a  received 
signal. 


With  respect  a  received  pulse,  we  assume  that  the  adaptive-thresholding  system  is 
capable  of  rejecting  any  level  crossing  which  occurs  too  prematurely,  that  is,  before  a  cer¬ 
tain  starting  time.  Under  this  assumption,  we  may,  without  loss  of  generality,  focus  atten¬ 
tion  on  a  single  received  pulse.  We  seek  to  determine  the  time  of  die  first  threshold  cross¬ 
ing  to  occur  after  a  starting  time  defined  to  bo  ;*  0.  A  practical  system  to  accomplish 
this  rejection  is  shown  in  Fig.  4.  The  combined  signal  and  noise,  r(f),  is  applied  to  a  fixed 


Fifj.  3—ExaRgerated  example  of  a  multiple  ihw>  hold  crossing 
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Fig.  4— System  for  the  elimination  of  false  alarms 


thresholder,  the  level  of  which  is  set  to  a  fixed  fraction  of  the  minimum  expected  amplitude 
of  the  received  pulses.  The  output  of  the  fixed  thresnolder  initiates  a  counter.  This  counter 
continues  to  count  as  long  as  the  Sxed  threshold  level  is  exceeded  by  r(t).  If  the  count 
reaches  a  certain  value,  the  gate  Is  enabled,  and  r(f)  is  allowed  to  pass  to  the  adaptive 
thresholder.  A  reasonable  critical  value  of  the  count  is  one-half  of  the  minimum  signal 
duration.  !f  the  minimum  signal  duration  is  large  enough  such  that  there  is  negligible 
probability  of  a  noise  pulse  exceeding  this  duration,  the  system  of  Fig.  4  will  reject  spurious 
pulses  occurring  outside  the  immediate  vicinity  of  the  signal  edge.  A  reset  pulse  prepares 
the  counter  for  the  next  received  signal.  In  keeping  with  our  strategy  to  detect  the  first 
level  crossing,  the  delay  in  the  upper  branch  should  exceed  somewhat  the  critical  value  of 
the  count. 


With  tc  a  random  variable  equal  to  the  first  threshold  crossing  time  of  a  received  pulse, 
4(f)  ib  defined  to  be  its  probability  density  function.  To  derive  an  expression  for  the 
probability  density,  we  relate  it  to  another  function  more  readily  obtained.  Referring  to 
Fig.  6,  let  P(t)e  represent  the  probability  that  a  first  threshold  crossing  of  r(f)  occurs  in. 
the  interval  between  f  and  f  +  e,  assuming  no  crossing  has  occurred  before  time  f.  From 
the  above  definitions,  it  fo&'ovs  (7)  that 


m  = 


4(f) 


Since  fc{t)  is  a  probability  density  function  and  we  have  assumed  4(0  =  0  for  f  <  0, 
it  is  required  that 


03 

fc{x)dx  =  1.  (2) 

Using  Eqs.  (2)  as  a  bcundry  condition,  Eqs.  (If  may  be  inverted  to  yield 

4(f)  ■  0(0  exp 


(3) 
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Fig.  5— Small  time  interval  of  a  received  pulse 


Once  an  expression  for  p(t)  is  deriyed,  Eqs.  (3)  will  provide  all  the  information  needed  for 
our  evaluation  of  the  performance  of  adaptive  thresholding. 

If  the  time  interval  between  t  and  t  +  e  is  chosen  small  enough  and  r(f)  is  assumed 
to  be  differentiable,  we  can  write 


r(f  +  e)  =  r(t)  +  r'(t)e,  (4) 

where  r'(t)  is  the  time  derivative  of  r(f).  This  equation  implies  that  we  may  ignore  the 
possibility  of  more  than  one  threshold  crossirg  in  the  small  time  interval.  We  now  assume 
that  r(t)  is  a  Markov  process.  In  the  definition  f  B(t)e,  it  follows  from  the  Markov  prop¬ 
erty  that  the  condition  that  no  crossing  occurred  beto*«-  t:me  t  is  equivalent  to  the  condi¬ 
tion  that  r(t)  is  less  than  the  threshold  at  time  f.  Thus  if  P[  ]  signifies  the  probability 
of  the  event  in  the  brackets, 

(S(t)e  =  P[r(t  +  e)  >  A\r(t)  <  A] 

-  P[r(t)  +  r'(t)e  >A>  r(t) ] 

P[r(t)  <  4] 

where  Eq.  (4)  has  been  used  in  the  last  step.  We  define  frr'(x,  y)  to  be  the  joint  probabil¬ 
ity  density  of  r(t)  and  r'(t).  The  probability  density  of  r(t)  is  represented  by  fr(x).  In  gen 
eral  the  densities  will  be  functions  of  time.  In  terms  of  the  density  functions,  Eqs.  (5) 
can  be  rewritten  as 


f  f  frr'(*.  y)  dxdy 

0  'A -ye 


P(t)e  = 


(6) 
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where  the  region  of  integration  in  the  nun:  .rator  can  be  determined  by  inspection  of  the 
rr  plane.  The  time  dependence  of  the  density  functions  is  implicit.  In  the  limit  of  small 
e,  the  inner  integral  of  Eq.  (6)  can  be  evaluated.  Taking  the  limit  of  Eq.  (6),  we  obtain 

f  yfrr'(A<  y)  dy 

■'o 

m  -  — - - .  (?) 

f  fAx)  dx 


s 

The  substitution  of  Eq.  (7)  into  Eq.  (3)  gives  a  complete  solution  to  our  problem.  How¬ 
ever,  the  result  will  seldom  be  a  simple  closed-form  expression,  since  the  density  functions 
of  Eq.  (7)  are  time  dependent,  making  the  evaluation  of  the  integral  in  Eq.  (3)  difficult 
in  most  cases. 

We  could  have  attempted  to  derive  fc(t)  from  the  Fokker-Planck  equation  for  con¬ 
tinuous  Markov  processes  and  first  passage  time  theory.  However,  the  present  approach  is 
easier  and  yields  greater  physical  insight 


STATIONARY  NOISE 


In  most  applications  the  noise  can  be  considered  stationary  over  the  observation  inter¬ 
val  In  this  case  the  general  results  can  be  simplified  somewhat.  Consider  the  received 
pulse  shown  in  Fig.  6.  Time  is  divided  into  four  periods.  Region  I  is  defined  by  0  <  f 
<<0>  the  time  before  the  occurrence  of  the  leading  edge.  Region  II  is  defined  as  the  inter¬ 
val  of  the  leading  edge,  t0  <  t  <  tv  Region  III  is  the  interval  between  the  leading  and 
trailing  edges,  <  t  <  t2.  Region  IV  is  the  remaining  time  period.  These  definitions  are 
not  precise  but  are  often  useful  when  the  pulse  has  a  reasonably  fast  rise. 


In  Region  I,  r(£)  has  no  signal  component,  by  definition.  Since  the  noise  is  station¬ 
ary,  neither  frr'(A,  y)  nor  fr(x)  has  a  time  dependence.  Thus 


m  »  c0. 


0  <t<t0. 


(8) 


and 


fc(t.)  o  C0<Tc°f  0  <t<t0,  (9) 

where  C0  is  a  constant.  Using  Eq.  (8)  in  Eq.  (3),  it  follows  that  for  Region  II 

,.t 


fJt)  =  e~c^m  exp 


j  m 

'to 


dt 


tn  ^  t  ^  1 1  . 


(10) 


in  Region  III  the  signal  component  of  r(t)  is  approximately  constant,  by  definition.  Once 
again  the  stationarity  of  the  noise  leaves  no  time  dependence  in  P(t).  It  is  then  easy  to 
see  that  fc(t)  has  the  form  given  by 

fc(t)  =  ClC'c2'. 


1 1  <  t  <  t2, 


(11) 
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where  C1  and  C2  are  constants.  In  most  applications  the  value  of  fc(t)  is  negligibly  swap 
in  Region  IV. 

The  general  features  of  the  variation  of  fc(t)  are  seen  in  the  previous  expressions.  We 
will  now  consider  some  important  special  cases. 


STATIONARY,  ADDITIVE  NOISE 

If  the  noise  is  additive,  we  can  write  r(f)  =  s(t)  +  n(t),  where  s(t)  is  the  signal  and 
n(t)  is  the  noise.  In  the  presence  of  additive  noise,  Eq.  (7)  becomes 


/, 


[y  +s{t)]fnn'{A  -s{t),'-)dy 


m  = 


-*(0 


I 


A-s(i) 


(12) 


fn(x)dx 


where  s'{t)  is  the  time  derivative  of  s(t),  fnn'(x,  y)  is  me  jv,r.t  probability  density  of  the 
noise  and  its  derivative,  and  fn(x)  is  the  probability  density  of  the  noise.  Note  that,  in 
contrast  to  Eq.  (7),  the  time  dependence  is  exhibited  explicitly  in  Eq.  (12).  Because  theQ,0 
noise  is  stationary,  there  is  no  implicit  time  variation  in  fnn'(x,  y)  or  fn(x). 

Equation  (12)  can  easily  be  evaluated  in  two  special  cases.  In  the  first  one,  it  is  as¬ 
sumed  that  P[n'(t)  <  -  s'(t )]  is  negligible  for  all  s(t)  and  n(t)  in  Region  II.  Formally,  we 
can  write 

fnn'(x,  V)  **  0,  y  <  -s'(0,  f0  <  *  <  u  •  (13> 

It  is  further  assumed  that  n'  is  symmetrically  distributed  about  the  origin,  that  is, 


fnn'(x,>y)  fnrt''X*  ~y)' 


(14) 
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Under  these  two  assumptions.  Eq.  (12)  is  approximately  given  by 

z*!u. 

rAs(t) 

I  fn(x)dx 

v,~oo 

Looking  at  the  right  side  of  Eq.  (15),  we  notice  feat  we  can  write 


(15) 


0(0-“  In 
dt 


~As{i) 


j  f„(x)dx 


(16) 


Using  the  £q.  (16),  j3(f)  can  be  readily  integrated.  Performing  fee  integration,  noting  that 
s(f0)  =  0  by  definition,  and  substituting  into  Eq.  (10)  yields 

e~e0fo  n 


f.(t)  =  s'(t)fn(A-8(t )) 


ft- 


(*)  dx 


t0*t<t  i, 


(17) 


where  the  factor  in  brackets  is  a  constant.  Equations  (9),  (11),  and  (17)  constitute  a 
close:! -form  solution  for  the  probability  density  in  this  special  case. 

For  the  second  special  case,  we  assume  s'(t)  is  a  constant  M  (not  necessarily  large) 
in  Region  II.  We  furtner  assume  that  n(t)  and  n'(t)  are  statistically  independent.  The 
latter  assumption  can  be  stated  formally  as 


fnn'&y)  =  fn(x')fn’(y) 

Under  these  two  assumptions,  if  M  4  0,  Eq.  (12)  reduces  to 

Mfn(A~s(t)) 


m  =  K 


f 


A-S(t) 


(18) 


(19) 


fn(x)dx 


where  If  is  a  constant  defined  by 


f 


-M 


(y)  dy,  M  4  0. 


(20) 


Clearly,  Eq.  (x9)  can  be  put  ir.  the  form  of  Eq.  (36).  It  then  follows  from  Eq.  (10)  that 


fed)  =  fn(A-s(t)) 


rA-S(t) 

K- 1 

KMe~CQt° 

1  fn(x)dx 

rA  \K 

mi 

fn(x)dx\ 

LV-~  /  J 

,  t0<t<t  i, 


(21) 
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where  the  last  factor  in  brackets  is  a  constant.  Equations  (9).  (11),  and  (21)  constitute  a 
closed-form  solution  for  the  probability  density  in  this  specie  case. 

It  is  noticed  that  in  the  limit  of  large  M,  Eq,  (20)  indicates  that  K  =  1.  Then  it  is 
seen  that  Eq.  (21)  reduces  to  Eq.  (17). 

It  can  easily  be  shown  (1)  that  additive  Gaussian  now  fulfills  the  conditions  specified 
by  Eq.  (14)  and  (18).  Thu?  if  Eq.  (13)  is  satisfied,  we  can  use  Eq.  (17)  for  the  Gaussian- 
noise  problem.  Alternatively,  if  M  is  constant,  we  can  use  Eq.  (21)  for  the  Gaussian-noise 
problem.  A  stationary  Gaussian  Markov  process  must  have  a  power  spectrum  of  the 
form  (7) 


Si  u) 


a 

CO  2  +  b  2 


where  a  and  b  are  constants.  A  flat  noise  spectrum  may  be  closely  approximated  over  any 
finite  frequency  range  by  a  suitable  choice  of  values  for  die  two  constants.  Of  course, 
stationary  white  Gaussian  noise,  which  is  Cat  over  the  infinite  frequency  range,  is  neither 
differentiable  nor  Markov. 

In  many  practical  situations  the  received  pulse  hrs  a  very  short  rise  time.  Rather 
than  use  the  formulas  already  established,  we  seek  a  direct  derivation  of  the  formulas  for 
this  special  condition 


DIRECT  DERIVATION  FOR  LARGE-SLOPE  LEADING  EDGES 

If  the  ieading  edge  of  the  received  pulse  has  a  large  slope,  it  is  intuitively  reasonable 
to  expect  that  there  is  only  a  small  probability  of  more  than  one  threshold  crossing  in 
Region  II  of  Fig.  6.  Formally,  if  fr'(y)  is  the  probability  density  of  the  slope  of  r(t),  we 
assume  that 


fr'iy )  «  0,  y<0,  (22) 

In  other  words,  the  signal  slope  is  sufficiently  great  it .  the  combined  signal  and  noise  to 
nave  negligible  probability  of  nonpositive  slope.  To  simplify  matters  further,  we  set  f0  - 
0,  that  is,  we  eliminate  Region  I. 

Let  Fr(x,  t)  be  the  probability  that  r(t)  <  x.  Thus  Fr(x,  t)  is  a  time-varying  probabil¬ 
ity  distribution  function.  Since  there  can  be  only  one  crossing  between  0  <  t  <  tx>  the 
probability  that  no  crossing  occurs  before  time  t  is  equal  to  Fr(A,  t),  for  if  r(f)  <  A , 
then  no  crossing  could  have  occurred  before  time  f.  Let  Fc{t)  represent  the  probability 
distribution  function  of  the  first  crossing,  that  is,  Fc(t )  is  equal  to  the  probability  that  a 
crossing  occurred  before  time  t.  From  this  definition  and  the  preceding  statements, 

1  -  Fc(t)  =  Fr(A,  t),  0  (23) 

Differentiation  of  a  distribution  function  givv.s  a  density  function.  Thus  at  all  points  where 
a  derivative  is  defined,  we  have  the  following  probability  density  function  for  a  threshold 
crossing: 
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fe{t)«--!j-Ft{A.t)  (24) 

We  may  ignore  the  possibility  of  a  first  crossing  for  t  >  /  j  if 

C‘  i 

1  -  fc(t)dt  «  1.  (25) 

JQ 

In  deriving  Eq.  (24),  we  have  not  invoked  the  Markov  property.  Thus  Eq,  (24)  is 
valid  for  processes  with  short  rise  times  and  differentiable  distribution  functions.  We  now 
turn  to  examples  of  the  applications  of  this  equation. 


10  DCN  J.  TORRJER! 

Equation  {30)  gives  the  condition  under  which  v/e  may  ignore  the  region  t  >  t1. 

We  can  easily  calculate  the  expected  value  of  lc,  the  random  variable  specifying  the 
pulse  position.  We  have,  if  Eq.  (25)  is  valid  and  /I  =  Mtv% 


E[U 


Using  Eq.  (29)  in  Eq.  (32),  we  obtain 

EUc)  23 


l 


r 


2  A/M 


tfc(t)  dt. 


(32) 


A  A 

—  1  -  2  erfc  i — 

M !  \a 

L  ' 


(33) 


Thus  the  arrival-time  estimation  by  adaptive  thresholding  at  the  half-amplitude  level  is 
biased.  However,  Eq,  (30)  indicates  that  the  amount  of  bias  is  small  and  that  this  bias  de¬ 
creases  n.pidly  as  signal  amplitude  increases. 

A  measure  of  the  fluctuation  in  pulse  position  is  provided  by  the  standard  deviation 
of  the  airivai  time,  which  is  defined  to  be 


s  -  v'eIP: 3  -  {SftTP  . 

When  Eq.  (25)  ic  satisfied,  we  can  write 

r2  A/M 

E[t* I  -  j 

VG 

From  Eqs.  (29)  and  (33)  through  (35),  it  follows  that  for  additive  Gaussian  noise, 


(34) 


(35) 


a\  // 

r 

i\  j  (a\ 

hg  2  erfc  ( - 

\Ml  \t 

U  j  W 

u  J 

a\~  -  M\  /ZlA\  -A 2\] 

—  1  -  2  srtc  ( — !  -  v' "  ~  wpl-r-s  !■ 

\M  [  \oJ  -*  ir\c  *  \2c‘/| 


1 36) 


Additive  Gaussian  noise  appears  in  the  models  for  coherent  amplitude-modulation 
(AM)  systems.  However,  noncoherent  AM  systems  usually,  contain  envelope  detectors.  At 
the  output  of  the  detector,  r{t)  is  described  by  a  Rice  distribution.  It  can  be  phown  that  the 
probability  distribution  function  js  (8) 

FM,  „  .  exp  «P  i-iU  (S’) 

\  2®  ./.A,  X  0  /  \  */ 

where  IQ{x)  is  the  zero-order  Bessel  function  of  imaginary  type.  From  Eqs.  (24)  and  (28), 
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where  / j(x)  is  the  first-order  Bessel  function.  The  substitution  of  Eq.  (38)  into  Eqs.  (32) 
and  (35)  results  in  integrals  which  cannot  be  evaluated  analytically.  However,  with  the 
aid  of  a  digital  computer,  it  is  possible  to  numerically  compute  the  value  of  tn,  the  nor¬ 
malized  expected  arrival  time,  defined  by 

(39) 

It  is  also  possible  to  determine  the  function  g(A/u),  which  is  defined  by 


M 

2a 


E[tci. 


(40) 


The  results  of  the  numerical  computations  are  summarized  in  Table  1.  The  first  column 
is  the  signal-to-noise  ratio  (S/N)  as  defined  in  the  usual  AM  system.  Thus 

Ez  2 A2 

S/N  -  lOlogiQ^  =  10  log  jo  -"pr*  (41) 

Hie  fifth  column  is  calculated  from  Eqs.  (36)  and  (40).  It  is  noticed  that  the  standard 
deviation  for  the  Gaussian  distribution  (coherent  AM  system)  is  uniformly  less  than  the 
standard  deviation  for  the  Rice  distribution  (noncoherent  AM  system).  For  a  fixed  value 
of  S/N,  the  degradation  in  standard  deviation  for  the  Rice  distribution  is  defined  tn  be 
the  additional  power  required  to  achieve  the  standard  deviation  that  exists  in  the  Gaussian 
case  for  the  specified  S/N.  The  sixth  column  lists  the  degradation  in  decibels  as  a  function 


Table  1 


3igna!~to-3cbe  Dependency  cf  Various  Quantities 


S/N 

(dfi) 

AJn 

m 

(Gauss) 

r~  1  . 

Approximate  j 
Begradftticm 
(db) 

12 

2.815 

0.462 

1.048 

0.987 

0.44 

14 

3.544 

0.478 

1.000 

0.24 

16 

4.462 

0.487 

1.015 

1.000 

0.13 

18 

5.617 

0.492 

1.009 

1.000 

0.08 

20 

7.071 

0.495 

1.005 

1.000 
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of  S/N,  For  S/N  <  12  d.B,  it  is  found  tlist  Eg.  (25)  is  not  well  satisfied  for  either  the 
Rice  or  the  Gaussian  distribution;  hence,  these  values  of  S/N  have  not  been  Fried. 

Examination  of  the  third  column  in  Table  1  shows  that  the  estimation  procedure  is 
biased.  However,  the  bias  steadily  decreases  as  S/N  increases. 


ASYMPTOTIC  RESULTS 

We  conclude  with  a  heuristic  discussion  of  the  asymptotic  properties  of  arrival-time 
estimation  by  adaptive  thresholding.  It  ha-  been  seen  in  the  last  section  that  the  esti¬ 
mation  procedure  is  biased.  Of  course,  some  of  the  apparent  bias  is  due  to  our  having 
assumed  a  negligible  probability  of  a  first  threshold  crossing  for  t  >  tx.  Some  bias  is  also 
due  to  setting  the  threshold  at  the  half-amplitude  level.  For  both  the  Rice  and  Gaussian 
distributions,  the  arrival-time  estimate  approaches  the  unbiased  value  of  one-half  the  rise 
time,  in  the  limit  of  large  S/N'.  Thus  we  can  say  that  the  estimate  is  asymptotically  un¬ 
biased. 

For  the  unbiased  estimate  of  arrival  time  of  a  low-pass  pulse  with  large  S/N  in  white 
Gaussian  noise,  the  Cramer-Rao  bound  requires  (6)  that 


>2  >  j*!  ; W'p 


(42) 


where  N  is  the  noise  power  spectral  density,  if  s'(0  =  M  is  a  constant  and.  tl  -  2A/M, 
Eq.  (42)  reduces  to 


s2  ^ 


N 

MU 


M2  ’ 


(43) 


where  we  have  assumed  that  the  rise  time  is  approximately  equal  to  the  inverse  of  the 
system  noise  bandwidth  B.  Thus  we  have  set 


N 

h 


NB 


(44) 


where  o2  is  the  mean  noise  power.  A  comparison  of  the  fifth  column  of  Teble  1  and 
Bq.  (40)  with  Eq,  (43)  indicates  that  the  adaptive-thresholding  estimate  approaches  the 
minimum  variance  expressed  by  the  Cramer-Rao  bound.  Thus  adaptive  thresholding  is 
asymptotically  efficient  for  white  Gaussian  noise  if  £q,  (44)  is  valid.  Clearly  the  perform¬ 
ance  of  a  noncoherent  system  is  bounded  by  the  performance  of  the  corresponding  co¬ 
herent  system.  Thus  from  a  comparison  of  tire  fourth  and  fifth  columns  of  Table  1,  it 
appears  that  adaptive  thresholding  is  also  asymptotically  e*.“‘cient  for  the  Rician  process. 
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